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Abstract 

We have looked at the evaluation of the riemann zeta function at odd 
arguments and have provided a simple formula to approximate the 
value with exponential convergence. We have compared it with var- 
ious other formulae present in literature. We have also evaluated an 
expression for the zeta function on the plane Sft(s) = 1. 



1 Introduction 

The Riemann Zeta function is a widely talked about problem which features 
as one of the Seven Millennium Problems. Simply Stated Riemann conjec- 
tured that 'All non-trivial zeros of the Zeta function have real part 
equal to one-half. If true, the hypothesis would have profound conse- 
quences on the distribution of primes in the integers |P88j . The Analytic 
Continuation of Dirichlet Series or the defining Riemann Zeta function is 
the simply given by 



CM 




K(s)>l 
U(s)>l,s ^ 1 



Its evident that C(s) is an infinite series which converges when the real part of 
s is greater than 1. It can indeed be extended meromorphically to the whole 
complex s-plane except for a simple pole at s = 1 with residue l(|T86j)Thc 
Riemann zeta function satisfies the functional equation 

COO = 2'tt- 1 sin (^) T(l - s) c(i - s), (1) 
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neering student at the National Institute of Technology (NIT), Karnataka(India). 
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The general formulae well known related to the Riemann Zeta function 

are 

C(2n, = (-ir+'^f^ (2) 

where Bm is a Bernoulli number; hence the values of zeta function vanishes 
at even negative integers because B m = for all odd m other than 1. No 
such simple expression is known for odd positive integers. 

The practical importance of this zeta function is in the distribution of 
primes. When we look at the distribution of primes between 10* and 10 1+J 
which i£N. If numbers of prime numbers lying in these different intervals is 
counted this distribution is related pretty closely to the zeros of the Riemann 
Zeta function. The distribution of primes is of prime concern in areas of 
security, encryption, factoring, etc. For e.g. the important RSA encryption 
scheme relies on the difficulty of detecting prime numbers. 



2 Uniform Convergence 

Uniform convergence is, for any e>0, there exists N such that for the same 
N, no matter which x you pick n>N implies \f n {x) — f{x)\<e 

The fundamental difference between point wise and uniform convergence 
can be best visualized by seeing that uniform convergence essentially implies 
that for large n, f n (x) is very close to f(x) as a graph everywhere, whereas 
point wise convergence simply says that tracking the value of a single point 
over time will give the right value. 

One example of a sequence that doesn't converge uniformly is x n on 
[0,1]. You can see this as, if x<l, x n converges point wise to 0, i.e. for a 
given point, it converges to zero. If x = 1, x n converges to 1 obviously. But 
for a given e, say 1/4, for any N we can find a point {\) l ^ N so if x is that, 
x N =l/2. 

The 3 known methods of proving uniform convergence are Abel's Uni- 
form Convergence Test, Weierstrass M-Test, Uniform Norm Method. We 
shall be employing the Uniform Norm Method to prove the uniform conver- 
gence of our functions. 



2.1 Uniform Norm Method 

If A C R, f:A — > R , / is bounded. Uniform Norm on A is 
||/|U : sup{\f(x)\:xeA} 

Theorem 1 A sequence (f n ) of bounded functions on A C R converges 
uniformly on A to f if and only if \ \f n — f\\^ — > 
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3 Zeta Function at Odd Arguments 



3.1 Intuition 

In recent years there has been a renewed interest in representing (,(2n + 1) 
(n€ N), N being the set of positive integers, by means of series which con- 
verge more rapidly than the defining series in Section [TJ These developments 
seem to have stemmed from the use of the familiar series representation for 
the known constants like C(3), C(5), C(7) in |SCjA00| . [VT0] 



C(3 ) = m (4) 

^ ' 7 ^ (2k + l)(2k + 2)2 2k y ' 

. 12 A 1 39^ 1 l^ l 

C v > 1^ n 5 s i nri (vrn) 20 ^ n 5 (e 2nn - 1) 20 ^ n 5 (e 2nn + 1) 

n=l n=l n=l ' 

19 °° 1 

C(7) = 06700^ _2 ^n7( e 2-_l) (5) 

n=l 

The main aim in the result presented below is to reduce the complexity of the 
huge computational formulae presented in literature. The intuition was that 
something simpler would exist which could at smaller values of argument 
might not be exact but with increasing values give better approximations 
to + 1)- An explicit mention has been made in the Section 13.31 with 
the other known results for computation of C,{2k + 1). In the |V10] Vestas 
has explicitly mentioned the computational difficulties of higher values of 
Bernoulli numbers. Thus we have also used an earlier result |A10] to avoid 
the need of Bernoulli numbers due to the computational factors as mentioned 
in [FP07] . 

3.2 Main Result 

In general it is known that 

, ^ 1 1111 

C(*)=>.Z7 = + ^ + + 



n=l 



1 1 1 1 1 1 1 1 1 

( ) + CgJ + 4I + 6*J + (3I + g7 + 27*^ + '5* + 25 1 + 125 1 

^ C(s) 1 1 

(1) + H 1 h • • • 

w 2 s 3 s - 1 5 s - 1 

1 + ^ + 1 



2 s k s - 1 
V l- 1 



Note: We would like to highlight that in the above formulation terms 
are fundamentally the poweres of primes and notions that terms such as 
15 or 21 which is a product of 2 primes would not repeat. Euler Product 
formula connecting the Primes and the Zeta function is 

P s 

C(^) — ^P>2,p is "prime~ s J" (10) 

Claim 1: Assuming^ is a set of primes, t(s)= J2 P >2 F=i > C(«)p>2 = 
and given the relation 



t(2s) _ t(2s + l) 

C(2s) J C(2s + l) 1 ' 

We can give a bound on the value of f(s) and present a convergence of the 
form 

lim \f(s) - 21 < e V e>0 (12) 

□ 

The above is the main crux and its pure observational that ' f(s) — > 2' . 
From Matlab if can be realized that the constant / which links both the 
relations averages with higher values of s is / = 2 . Hence, initially when s 
is small it is different from 2 but as s increases to above 20 the value of / 
averages to 2. Hence during computation of £(s) for higher values of s we 
can use the value / = 2 and use the formula. 



t(2s + l) . . . . 

C 2s + 1 =2 \ (9 ; > C(2s) 13 
t{2s) 

Using Claim 1 with the notations of t(s) and and £(s) remaining consis- 
tent. 

E^T = C(*)(l-2- s )-l (14) 

p>3 P 

E^ T (=^))-^ L T = C00(l-2-*)-l (15) 



p>2 ^ 



t(a) = C(^)(1-2" S )-1 + ^ T (16) 
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Using the above we can rewrite Equation [IT] 

t(2s) _ C(2 S )(i-2- 2s )-i+ i2 i- T 



((2s) 



C(2 S ) C(2s)(2^-1) 

t(2a + l) c(2 S +i)(i-2-( 2s + 1 ))-i+ 525 ^ rzT 



(17) 



(19) 



((2s + 1) " « 2s+1 ) 

= 1 - 2~( 2s+1 ) - <^5+T) + C(2s+l)(W i -l) ( 20 ) 
Using Claim 1 we can rewrite the above as 

t(2s) t(2s + l) 

((2s) ~ ; C(2s + l) 
(l _ 2 " 2 * - _L- + __I ) = / / 1 _ 2 -(2,+i)_ 1 + 



C(2s) C(2s)(2 2 * - 1)/ J \ " C(2s + 1) C(2s + l)(2 2s+1 -l)- 



2 — 2 2 s 

c(2*+i) = I ;_ 2S 1 r ( 21 ) 



(2.2- - 1)( 2a ^M"- 1 -1 + 5^ 

_ (-l)- +1 B 2s (2^) 2 - 
~ 2(2s)! we & et 

expression for C(2s + l)Vs > 1 



Using the formula for £(2s) = - ^ 2 (2V)1 — we the following final 



, f2s n /.2 2 ^(l-2 2 *) 

^ J fo2.s+1 4(2s)!(l-2^-i) 2 ^+i(l-/)+/-2 \ 1 J 

^ Z ^ V(-l) s+1 S 2s 7r 2s (2 2 '>-l) 2 y 

In order to avoid problems with huge Bernoulli results we could use the 
Result in |A10j and rewrite the Equation [22] as 

«*) = {2k _ ^ _ !„ <«>■ + 2) + (^)^(-D' 

C(2fc + 1) = -/- 2(1-^)2(2^ I ^ 



( 2 " 2fc ) " 22fc+1(1 2 /)+/ ' 2 ES^pd^C^ + 2) + ^ 



Numerical Results to authenticate the above formula (Using / = 2) || 

We have used the value of / = 2 for computing the table below. This provides us 
accurate results for higher values of computation of f (s) 
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Formula Result 


Actual Result 


Difference 




C(3) 


1.21992 


1.202056 


1.7861 x 10' 


-2 


C(5) 


1.03933 


1.036927 


2.3021 x 10' 


-3 


C(7) 


1.00861 


1.008349 


2.4187 x 10- 


-4 


C(9) 


1.00204 


1.002008 


2.5985 x 10" 


-5 


C(n) 


1.00494 


1.004941 


2.8476 x 10" 


-6 


C(13) 


1.00012 


1.000122 


3.1468 x 10" 


-7 


C(15) 


1.00003 


1.000030 


3.4890 x 10" 


-8 



3.3 Notes 

There are plenty formulae for Riemann zeta numbers at odd numbers but 
what differentiates the above from the rest of them are: 



• The convergence is Exponential in Nature. We can see that the dif- 
ference of the actual zeta values at odd numbers and the zeta values 
according to the formula proposed by us. We state that with higher 
values of zeta the difference would be almost negligible with orders of 

o(io- fe ). 

• The Simplicity in evaluation. A few formulae presented in literature 
from jZW93j . [SGAH0] . [S98] are as follows : 



C(2ro + 1) 



(-l) m 7T 



m _2m 



log(2) 



(1 - 2- 2m ) V T(2m + 2\ 



+E 



(2 



>l-2n 



)r(2n)C(2n) 



71=1 



r(2m + 2n + 2) 



+ 



1 - 2" 2m 



m— 1 

E 

n=l 



(2 



In— 2m 



1)- (ir 2 ) n ((2m - 2n + 1) 
T(2n + 2) 



n € N (23) 



+ < 2 ") ! Ep^(|^)cw + 1 >)»^ w 
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n-1 



(2»)l£ 



(3 2j - 1 



C(2j + 1) 



(2n) 

7(3) ^ (2n - 2j + 1) 



E 



.J x (2n-2j)! V (2vr)2i 

2C(2j,i)-(3^-l)C(2j)x 



(27T) 



2J-1 



J?iGN (25) 



C(2n + 1) = (-l)"-i 7 — - flog(2) 



(2n)!(2 4 ™+! + 2 211 - 1) 

+ E (^W + (2n) ! E T^fe &) ^ + x > 

(2r yy (ZIZ C(2^)-2^(2^-l)C(2 J ) V 
^•Z>( 2n _2j + 1)!- (2^-1 J nG ^ 

(26) 

Quite explicitly it can be seen that our formula is independent of the 
/mge number of double sums and integrals which occur in the other 
evaluations of C(2n + !)• The convergence of the formula mentioned 
might be exponential and little slower than other formulae present in 
literature but we feel reducing complexity(by reducing the parameters, 
increasing ease of calculation by softwares like Mathematica, Matlab, 
Maple etc.) helps in the same evaluation time of higher values of zeta 
functions at odd arguments. We also negate the need for Bernoulli 
Numbers which at higher values might be huge and cause extra evalu- 
ation time. Time complexity of Bernoulli Numbers has been explicitly 
mentioned in |FP07j . In a certain ways our expression also resembles 
a few of them presented above. 

Preknowledge of Odd £ Values isn't Required: It can be seen in the 
most of the previous formulae the sums include terms which require 
evaluation and hence pre-knowledge of £ at odd values. Thus our 
formula can be used to evaluate any C(k) without the need for knowing 
the values of £ (2 k — 1) for all values of k. 

Exact expressions for £ values at odd arguments are generally known 
for the famous ((3), C(5), C(7) where the convergence is tremendously 
fast and accurate. As such most of the formulae for general expres- 
sions of C(2& + 1) are mainly of approximations. And we would like 
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to mention that ours is one among the many approximations with 
exponential convergence and easier computation. 



• The presented formulation doesnt converge rapidly for higher values 
of zeta function purely because of the convergence of the of the zeta 
function to 1 with higher parameters, cause in that situation the con- 
vergence would uniform convergence and not exponential in nature. 
The exponential nature of the convergence suggests that the presented 
formula is an alternative expression for the zeta function at odd values. 

4 Riemann Zeta function at Re(s) = 1 

4.1 Intuition 

In general it is known that £(s) diverges for s = 1 and there are many 
definitions which meromorphically extend it to the entire plane except at 
the pole. We wanted to investigate the nature of the nature £(s) for values 
of the form s = 1+ib. It obviously isnt an unbounded value and turns out to 
be bounded. There are quite a few stringent uniform continuity conditions 
which are required to be tested for the below result, mentioned in Lemma f4.2l 
and 14.21 There has not been many known results exsistent for the £(1 + ib) 
and hence we thought it would help in looking at such a computation. 

4.2 Background Results 

— ib 

Lemma 1 f n {x) = x ^ x+n a^ is uniformly convergent except the singularity 
at {0} 

Proof of 1: Our proof shall be for a = 1 cause we are looking at £(1 + ib). 
The Uniform Norm Test in Section 12.11 cannot be applied because x varies 
from to oo in the integral and in the interval (0,1), - turns out to be 



unbounded. Thus we proceed as follows: f n (x) = ^L'+nl = ^ifj+n)^ ~ 



both a n (x) and b n (x) should be uniformly convergent. We shall divide the 
above into 2 intervals 

• x 6 (1, oo) 




a n (x) = 



cos (ft ln(x)) 
x(x+n) 



< 1 for all x, n > 




x(x+n) 



— x(x+n) 
1 _ r 



Vx > 1. Hence it is bounded. 



Since lim 



■71— >oo SUp 
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Thus 



and | \a n — a\ 







cos (ft ln(x)) 
x(x+n) 

Similarly the same is true for b n (x) , hence / n (x) is uniformly conver- 
gent in the interval (l,oo) 



i£ (0, 1) we rewrite the 

a n (x) 



cos(ft ln(a;)) 
x(x+n) 

Er = o(i-^) fc ^(foi"W) 



(27) 
(28) 



Thus as limrwoo, the singularity has been removed and the (} cos ( fel °W) 
converges uniformly in the compact subset of any plane not containing 
0. The above follows for b n {x) and thus f n (x) is uniformly convergent 
in the range (0,1] 



Hence the uniform convergence of the J ( 



oo x - 



x(x+n) 



gives us the option of ex- 



changing the summand and the integrand in Equation 14.31 □ 



Lemma 2 Both the functions (i)f\ n (x) = I ^ — — 



\^ 2n+x+l 

are uniformly convergent in the entire plane. 



, (ti)hn(z) — i(2n+x+l)(2n+x+2) 



Proof of 2: (i) f ln (x) = ( 

2n+z+i) fc — ^ Hence fi n (x) is bounded, 
lim^oo f ln (x) = = / V k G [2, oo) x G IR + + {0} and n G N + {0} 



\\fm(x)-f\\ 



-1 



2n + x + 1 



sup 

( — )* 
V 2ra 



1 



R+ 
A' 



2n + x + 1 



: x > 



We can infer that lim n ,_ >00 (^-) 2 ' c VA; G [2, oo) 



Hence 



fin 



/|| —7-0 . Thus /i n ( x ) is uniformly convergent to 0. 



(29) 
(30) 
(31) 



(H)f2n{x) 



X 



-ib 



(2n + x + l)(2n + x + 2) 

Q—ib ln(x) 

(2n + x + l)(2n + x + 2) 

cos(61n(x)) 
(2n + x + l)(2n + x + 2) 

cos(ft ln(x)) 



sin(61n(x)) 



(2n + x + l)(2n + x + 2)- 



a n(%) = { (2n+x+l)(2n+x+2) ) <h HeI1Ce f*>( x ) is bounded. 

limn-**, / 2n (x) = = / Vx G R+ + {0} and n G N + {0} 
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\a n (x) - a|| R 4 



cos(6 ln(x)) 



sup< 



(2n+x+l){2n+x+2) 
cos(b ln(a;)) 



(2n+x+l){2n+x+2) 



X > 



cos (ft ln(x)) 



< 



(2n+l)(2n+2) — (2ra+l)(2n+2) 



(32) 
(33) 
(34) 



Since lim r 



i 



Thus | \a n {x) 



0. The same follows for 



>DO (2n+l)(2n+2) 

b n (x). Hence f n (%) is uniformly convergent in the entire plane. □ 



4.3 Main Result 

The general expression for £(s) continued Re(s) > would be 



i 00 (— 



n=l 

Let s=a+ib Va>0,i^0 



(_-\\n rl — tb 

n=l 

Our aim or trick in handling the above expression is on the lines of [TH07| 
where we would like to bring an integral within the summation. We would 
like to highlight that the following integral would help our expression 35. 

x + n a dX ~ sinh(fg) n a ^ 

oo x -i -ib (n- 1 )- ib , . 

dx = ■ — TT 7 ! — \ " — 38 

x + n smn(zo7r) n 

Using the above in Equation 35, we get an expression as follows 

c(i + ^b) = - 1 , sinh( ^ } y ( -ir T (39) 

SV ' l-2~* b vr ^ Jo x + n y J 

Using Lemma 14.21 

roo — 1_,— ifc 
„=1 ^ x + n 

-1 sinh(i6vr) f°° / ^ (-l^x' 1 



c(i+») = - 1 .^vi-irf— ^ (40) 

71=1 U 

«r(f^)A (4i) 



1 - 2 



10 



We note that 

x + n 



n=l 



where the digamma function \l/(a;)= 

c(i + *) = : -! ,, sin y 67r) r n~ + 1) - n^v^x (43) 



1 - 2- lb 2ir 
We know the equation given in |GR07| 

oo 

tf(i + a)-tf(a) = 5Z(-l) fe C(^a)* fc_1 V|i| < |a| (44) 

C(k,a) = — (45) 

With a = 2±i and t=±. With x > . V|t| < \a\ Hence 

*(f + 1) - *(^) = E(-l) fc C(fc,^-)(^ 1 (46) 

k=2 

oo oo ^ 

~ ^ (2n + x + l)fc 

= yy 22 ^ (48) 



(2n + x + l) fc 

fc=2n=0 v ' 



^-^ In + x + 1 

k=2 n=0 



Using the above in Equation |3T] and Lemma 14.21 we get the following 



oo 00 00 , k 



JU fc=2n=0 

-1 sinh(ifeTr) f 00 A 
— lini / 7 dx 

l-2~ ib vr fc->oo J ^ (2n + x + l)(2n + x + 2) 

With Lemma 14.21 we can simplify the above to 

C(l + ib) = r :l lb SinH f 7r) f> + I)"* ( - incsch(inb) 

n=0 

+ mcsch{mb)((2n + l) i6 (2n + 2)" i6 )) 

(50) 
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-1 



oo 



C(l + »6) 



1 - 2~ ib 



^(-(2n+l)- l6 + (2n + 2)" i6 ) (51) 



n=0 



-1 



(-i) 



■it 



(52) 



1 - 2~ ib 



n=l 



Basically the above result stems from the first step but theres a subtle 
factor of n vanishes in the denominator which emphasis that the above 
is the actual formula for £(1 + ib). The rigorous need for proofs of uniform 
convergence show the innate siginificance of convergence when related to 
zeta(s). We could validate few known results like lim n _ >00 (s — l)C(s) = 1 > 
lim^o C(l + ib) — > oo , there is also the mention that the other zeros would 
lie on b = j^pj for any k ^ € I ( which can be got by simply relating the 
analytic continuation of Dirichelet Series and the Alternating Zeta function) . 
Hence we could find more identities with regard to the zeta function on the 
Line s = 1 as shown in [S03 

5 Conclusion 

This article gives a look into the evaluation of the Riemann Zeta func- 
tion with odd arguments and at the strip 3?(s) = 1 such that we could 
achieve more results and analyze the regions with better knowledge. The 
advantages of each formulation have been mentioned in each section. We 
might compromised a little on the accuracy but our main goal has been to 
reduce complexity(and providing exponential convergence) so that the out- 
come would be almost accurate but evaluated faster with lesser parameters 
and fewer huge calculations. We also have given a result for the £(1 + ib) so 
that better analysis could be made with the formula. 
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